A stochastic limit approach to the SAT problem 
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I. INTRODUCTION 

Although the ability of computer is highly progressed, 
there are several problems which may not be solved ef- 
fectively, namely, in polynomial time. Among such prob- 
lems, NP problem and NP complete problem are fun- 
damental. It is known that all NP complete problems 
are equivalent and an essential question is whether there 
exists an algorithm to solve an NP complete problem in 
polynomial time. They have been studied for decades 
and for which all known algorithms have an exponential 
running time in the length of the input so far. The stan- 
dard definition of P- and NP-problems is the following 
[11,15,14,18]: 

Definition 1 Let n be the size of input. 
( 1 )A P-problem is a problem whose time needed for solv- 
ing the problem is polynomial time of n. Equivalently, it 
is a problem which can be recognized in a polynomial time 
of n by deterministic Turing machine. 
(2) An NP-problem is a problem that can be solved in poly- 
nomial time by a nondeterministic Turing machine. 

This can be understood as follows: Let consider a prob- 
lem to find a solution of / (x) = 0. We can check in poly- 
nomial time of n whether xq is a solution of / (x) = 0, 
but we do not know whether we can find the solution of 
/ (x) — in polynomial time of n. 



Definition 2 An NP-complete problem is 
polynomialy transformed NP-problem. 



problem 



We take the SAT (satisfiable) problem, one of the NP- 
complete problems, to study whether there exists an al- 
gorithm showing NPC=P. Our aim in this paper and the 
previous papers [7,9,10] is to find a quantum algorithm 
solving the SAT problem in polynominal time of the size 
of the problem. 

Let X = {xi, • • • , x n } be a set. Then Xk and its nega- 
tion Xk(k — 1, 2, • • • , n) are called literals and the set of 
all such literals is denoted by X' = {x\, x\, • • • , x ni x n }. 
The set of all subsets of X' is denoted by TiX') and 
an element C G J- {X') is called a clause. We take a 
truth assignment to all Boolean variables Xk- If we can 
assign the truth value to at least one element of C, then 
C is called satisfiable. When C is satisfiable, the truth 
value t (C) of C is regarded as true, otherwise, that of C 
is false. Take the truth values as "true false <->0". 
Then Cis satisfiable iff t (C) = 1. 



Let L = {0, 1} be a Boolean lattice with usual join V 
and meet A, and t (x) be the truth value of a literal x 
in X. Then the truth value of a clause C is written as 
t(C) ee V xec t(x). 



Moreover the set C of all clauses Cj (j = 1,2, 



,m) 



is called satisfiable iff the meet of all truth values of Cj 
is 1; t{C) ee Af =1 t{C 3 ) = 1. Thus the SAT problem is 
written as follows: 

Definition 3 SAT Problem: Given a Boolean set X = 
{xi, ■ ■ ■ ,x n }and a setC = {C\, ■ ■ ■ ,C m } of clauses, deter- 
mine whether C is satisfiable or not. 

That is, this problem is to ask whether there exists 
a truth assignment to make C satisfiable. It is known 
in usual algorithm that it is polynomial time to check 
the satisfiability only when a specific truth assignment 
is given, but we can not determine the satisfiability in 
polynomial time when an assignment is not specified. 

Ohya and Masuda pointed out [7] that the SAT prob- 
lem, hence all other NP problems, can be solved in poly- 
nomial time by quantum computer if the superposition 
of two orthogonal vectors |0) and |1) can physically de- 
tected. This result was rewritten in [16] showing that 
OM SAT-algorithm is combinatoric. 

The output of the OM quantum-SAT algorithm is a 
superposition vector a |0) , and, in order to effec- 

tively implement this algorithm, it is necessary to dis- 
tinguish this superposition from the pure vector |0). If 
/3 is not zero but very small this detection is considered 
not to be possible with the present technology. In [10] it 
is shown that such a distinction can be realized by com- 
bining a nonlinear chaos amplifier with the OM quantum 
algorithm, which implies the existence of a mathematical 
algorithm solving NP=P. It is not known if the amplifica- 
tion method of Ohya and Volovich is in the framework of 
quantum Turing algorithm or not. So the next question 
is (1) whether there exists a physical realization combin- 
ing the quantum-SAT algorithm with chaos dynamics, 
or (2) whether there exists another method to achieve 
the above distinction of two vectors by a suitable uni- 
tary evolution so that all processes can be discussed by 
a quantum Turing machine (circuit). In this paper, we 
argue that the stochastic limit [1], can be used to find 
another method to solve problem (2) above. 

In Section 2 we review the mathematical frame of the 
OM quantum-SAT algorithm and in Section 3 we review 
the representation of this algorithm given by Accardi and 
Sabaddini [16]. In Section 4, we state the problem to 
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distinguish two vectors with a quick review of OV-chaos 
algorithm. In Section 5 the new notion of quantum adap- 
tive stochastic system is proposed and we show that it 
can be used to solve the problem NPC=P. The details 
(e.g., proofs) of this paper is discussed in [3]. 



II. QUANTUM ALGORITHM 

The quantum algorithms discussed so far are rather 
idealized because computation is represented by unitary 
operations. A unitary operation is rather difficult to re- 
alize in physical processes, more realistic operation is one 
allowing some dissipation like semigroup dynamics. How- 
ever such dissipative dynamics destroys the entanglement 
hence they very much reduce the ability of quantum com- 
putation of preserving the entanglement of states. In 
order to keep the high ability of quantum computation 
and good entanglement, it will be necessary to introduce 
some kind of amplification in the course of real phys- 
ical processes in physical devices, which will be similar 
to the amplication processes in quantum communication. 
In this section, to search for more realistic operations in 
quantum computer, the channel expression will be useful, 
at least, in the sense of mathematical scheme of quantum 
computation because the channel is not always unitary 
and represents many different types of dynamics. 

Let Ti be a Hilbert space describing input, computa- 
tion and output (result). As usual, the Hilbert space 
is U = ®f C 2 , and let the basis of H = ®f C 2 be: 



e (=|0» = |0>®-- 
|0)® |1> ■ ■ ,e 2 «-_i 
Any number t (0, 



® |0) ® |0),ei (= |1» = |0) ® ••• ® 
= (2^-1)) =|1)®---®|1>®|1>. 
■ • • , 2 N — l) can be expressed by 



N 



(fe) 
a t 

by 



or = 1, so that the associated vector is written 



|t)(=e t )=®£Li a[ k) 



And applying n-tuples of Hadamard matrix A = 
^ ^ to the vacuum vector |0) , we get 



Put 



A|0>(=£(0)) = ®f-^(|0> + |l>). 



Then we have 



J2 ex P ( ^Sr ) ' 

k=o x " 



which is called Discrete Fourier Transformation. The 
combination of the above operations gives a unitary op- 
erator U F (t) = W (t) A and the vector £ (t) = U F (t) |0) . 



A. Channel expression of conventional unitary 
algorithm 

All conventional unitary algorithms can be written as 
a combination of the following three steps: 

(1) Preparation of state: Take a state p (e.g., p = 
|0) (0|) and apply the unitary channel defined by the 
above Up (t) : A* F = Adu F (t) 



A* 



Ad Up => A* F p = U F pU* F 



(2) Computation: Let U a unitary operator on Ji rep- 
resenting the computation followed by a suitable pro- 
gramming of a certain problem, then the computation 
is described by a channel = Adjj (unitary channel). 
After the computation, the final state pf will be 



PS 



AuA 



u a fP- 



(3) Register and Measurement: For the registration of 
the computed result and its measurement we might need 
an additional system IC (e.g., register), so that the lifting 
£* n from S (H) to S (H ® JC) in the sense of [2] is useful 
to describe this stage. Thus the whole process is wrtten 
as 

Pf - £* m {A\jA F p) . 

Finally we measure the state in JC: For instance, let 
{Pk\k € J} be a projection valued measure (PVM) on 
K 



£(t) = W(t)£(0) = 



after which we can get a desired result by observations 
in finite times if the size of the set J is small. 



B. Channel expression of the general quantum 
algorithm 

When dissipation is involved the above three steps have 
to be generalized. Such a generalization can be expressed 
by means of suitable channel, not necessarily unitary. 

(1) Preparation of state: We may be use the same 
channel = Adjj F in this first step, but if the number 
of qubits N is large so that it will not be built physically, 
then A F should be modified, and let denote it by A* P . 
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(2) Computation: This stage is certainly modified to a 
channel Aq reflecting the physical device for computer. 

(3) Registering and Measurement: This stage will be 
remained as aobe. Thus the whole process is written as 

Ps = £* m (A* c A* pP ) ■ 



III. QUANTUM ALGORITHM OF SAT 

Let and 1 of the Boolean lattice L be denoted by the 

vectors |0) = ^ q an d |1) = ^ 1 ) * n ^ e Gilbert space 

C 2 , respectively. That is, the vector |0) corresponds to 
falseness and |1) does to truth. This section is based on 
[7,16,3]. 

As we explained in the previous section, an element x G 
X can be denoted by or 1, so by |0) or |1) . In order to 
describe a clause C with at most n length by a quantum 
state, we need the n-tuple tensor product Hilbcrt space 
H = ®"C 2 . For instance, in the case of n = 2, given 
C = {xi,x 2 } with an assignment x\ = and x 2 = 1, then 
the corresponding quantum state vector is |0) ® |1) , so 
that the quantum state vector describing C is generally 
written by |C) = \x\) <g> \x 2 ) G H with ij = or 1 
(k=l,2). ' 

The quantum computation is performed by a uni- 
tary gate constructed from several fundamental gates 
such as Not gate, Controlled-Not gate, Controllcd- 
Controlled Not gate [20,8]. Once X = {xi, ■ ■ ■ ,x n } and 
C = {Ci, C*2, • • • , C m } are given, the SAT is to find the 
vector 

\t(C))=Af =1 V xeCj t(x), 

where t(x) is |0) or |1) when x = or 1, respectively, and 
t(x) A t(y) = t(x A y), t(x) V t(y) = t(x V y). 

A. Logical negation 

Definition 4 Let X be a set. A negation on X is an 
involution without fixed points, i.e. a map x G X x' G 
X such that (x 1 )' = x ;x^x' Vi £ I, i' is called 
the negation of x. 

Proposition 5 Given a nonempty set X with a negation 
(x \— > x') and denoting, for / CI 

/' := {a/ G X : x G /} 

there exists a set I C X such that X = I U /'. 

Thus a finite set with a negation must be even. Let 
A be a finite set with 2n elements and with a negation 
(x i ► x'). A partition X = I U /' , |7| = n can be 
constructed with an n-step algorithm. 



Not all n-step algorithms are equivalent. The problem 
may come from the following question: having produced 
k elements of a given set, how difficult is to produce a 
mean element Xk+i, different from the previous ones? 

Definition 6 Given a set X with a negation x i— ► x' , a 
"clause" is a subset of X . A minimal clause is a subset 
I C X such that I (1 1' = <f> (i.e. if I contains x, it does 
not contain the negation of x). 

In a set X of cardinality 2n there are 2" minimal 
clauses. Given a set Co of clauses, if there are non min- 
imal clauses in it, then we can eliminate then from Co 
because any truth function must be identically zero on a 
non minimal clause. 

However, to eliminate the non minimal clauses from 
Co, one has to "read" all its elements. These can be of 
order 2™. 

B. Truth functions 

The set {0, 1} is a boolean algebra with the operations 

eye 1 :— max{e, e'}, e A e' :— min{e, e'} 

(e, e' G {0,1}). A clause truth function on the clauses 
on the set {x\, . . . , x n , x\, . . . , x' n } is a boolean algebra 
homomorphism 

t : Parts of {xi, . . . , x n , x\, . . . , x' n } — > {0, 1} 

with the property (principle of the excluded third ): 

t(x j )Vt(x' j ) = l; Vj = l,...,n (I) 

Because of (I), such a function is uniquely determined 
by the values {t(x\), . . . , t(x n )}, hence such functions are 
2". For this reason, in the following, we will simply say 
truth function on {x\, . . . , x n } meaning by this a truth 
function on the clauses of the set {x\, . . . , x n , x[, . . . , x' n }. 
Conversely given any n-ple e = (ei, . . . , e n ) € {0, 1}", 
there exists only one truth function on {x\, . . . , x n }, with 
the property that 

t( x j) = £ j 5 Vj = 1,.. .,n 

In the following, given a truth function t, we will denote 
e t the string in {t(x\), . . . , t(x n )} uniquely associated to 
that function. 

Let T be the set of truth functions on {xi, . . . , x n }. 
The function 

teT^ \t( Xl ),...,t(x n )) g ®"C 2 

defines a one-to-one correspondence between T and the 
set {0, 1}, that is, a one-to-one correspondence between 
truth functions and vectors of the computational basis of 
(g)"C 2 ) 
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Proposition 7 Let C C X be a clause and I, I' the sets 
associated to it through the procedure explained in § (1). 
Let t be a truth function on {x\, . . . ,x n }. Then 



t(C) 



Liei 



jei> 



Therefore as stated in Introduction, a set of clauses Co 
is said to be SAT if there exists a truth function t, on 
{xi, . . . , x n } such that 



t(C ) :=t{A C ec a C) 



n w 

cec 



1 



C. Quantum algorithm for the SAT Problem 

We review here a technique, developed in [7], which 
shows that the SAT problem can be solved in polynomial 
time by a quantum computer. 

Given a set of clauses Co = {C\, . . . , C m } on X, Ohya 
and Masuda construct a Hilbert space H = ®" +AI C 2 
where p is a number that can be chosen linear in mn, 
and a unitary operator U~c ■ H — > H with the property 
that, for any truth function t, 

Cco|et,0„ >= |et,a^_i,t(C ) > 

where, e t is the vector of the computational basis of 
<E> n C 2 corresponding to t, and 0^ (resp. x e x ) is a string 
of p zeros (resp. a string of (/i — 1) binary symbols de- 
pending on e). 

Furthermore Uq is a product of gates, namely of uni- 
tary operators that act at most on two q-bits each time. 

Let Co and Uc be as above and, for every e e {0, 1}™, 
let t e be the corresponding truth function. Applying the 
unitary operator Uc to the vector 



ee{0,l}" 



one obtains the final state vector \vf > 



\Vf >■= U Co \V) = ^ l £ ' *e( C 0) > 



Theorem 8 Co is SAT if and only if, 

P n+l >,iU Co \v) + 
where P n +^,i denotes the projector 

Pn+^,l ■= ln+^-l ® I 1 >< !| 

onto the subspace ofH spanned by the vectors 



According to the standard theory of quantum measure- 
ment, after a measurement of the event P n + M ,i, the state 
p = \vf >< Vf \ becomes 



Pn+fj,,lpPn+fi,l 

Trp'P n+ ^ 



=■ P 



Thus the solvability of the SAT problem is reduced to 
check that p' ^ 0. The difficulty is that the probability 



of P, 



Trp'P n+ ^ = \\P n+ ^f 



\T{C ) 
2" 



where |TCo)| is the cardinality of the set T(Co), of all the 
truth functions t such that t(Co) = 1. 

We put q := with r := \T(Co)\ in the sequel. Then 
if r is suitably large to detect it, then the SAT problem 
is solved in polynomial time. However, for small r, the 
probability is very small and this means we in fact don't 
get an information about the existence of the solution of 
the equation t(Co) — 1, so that in such a case we need 
further deliberation. 

Let us simplify our notations. After the quantum 
computation, the quantum computer will be in the state 



\vf) = Vl - I 2 \<Po) ® 1°) +Q\Vi) ® I 1 ) 

where \ipi) and \<po) are normalized n qubit states and 
q = yjr/2 n . Effectively our problem is reduced to the 
following 1 qubit problem. Wc have the state 



IV;) = v / l^|0>+9|l) 

and we want to distinguish between the cases q = and 
q > (small positive number). 

It is argued in [13] that quantum computer can speed 
up NP problems quadratically but not exponentially. 
The no-go theorem states that if the inner product of 
two quantum states is close to 1, then the probability 
that a measurement distinguishes which one of the two 
it is exponentially small. And one could claim that am- 
plification of this distinguishability is not possible. 

The proposal of [10] is that we do not make a mea- 
surement, which will be overwhelmingly likely to fail, but 
to use the output I of the quantum computer as an 
input for another device which uses classical chaotic dy- 
namics to amplify the required probability. 

Such a chaos amplifier introduces a new ingredient in 
quantum computation and in [9,10] it is proved that such 
an amplification is possible in polynomial time. 

In [9,10] a practical realization of the new amplifier 
mechanism is not suggested, however it seems to us that 
the quantum chaos amplifier considered in [10] deserves 
an investigation and has a potential to be realizable. 
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D. Chaotic dynamics 

Various aspects of classical and quantum chaos have 
been the subject of numerous studies, see [17] and ref's 
therein. Here we will argue that chaos can play a con- 
structive role in computations (see [9,10] for the details). 

Chaotic behavior in a classical system usually is con- 
sidered as an exponential sensitivity to initial conditions. 
It is this sensitivity we would like to use to distinguish 
between the cases q — and q > from the previous 
section. 

Consider the so called logistic map which is given by 
the equation 

x n+1 = ax n (l - x n ) = f(x), x n e[0, 1]. 

The properties of the map depend on the parameter a. 
If we take, for example, a — 3.71, then the Lyapunov 
exponent is positive, the trajectory is very sensitive to 
the initial value and one has the chaotic behavior [17]. 
It is important to notice that if the initial value x = 0, 
then x n = for all n. 

It is known [19] that any classical algorithm can be 
implemented on quantum computer. Our quantum chaos 
computer will be consisting from two blocks. One block 
is the ordinary quantum computer performing compu- 
tations with the output \ip) = \/l — q 2 |0) + <?|1). The 
second block is a computer performing computations of 
the classical logistic map. This two blocks should be 
connected in such a way that the state first be trans- 
formed into the density matrix of the form 

p = q 2 Pi + (1 - q 2 ) P 

where Pi and Pq are projectors to the state vectors |1) 
and |0) . This connection is in fact nontrivial and actually 
it should be considered as the third block. One has to 
notice that Pi and Po generate an Abelian algebra which 
can be considered as a classical system. In the second 
block the density matrix p above is interpreted as the 
initial data po, and we apply the logistic map as 

_ (I + f m (po)cr 3 ) 

Pm~ 2 

where / is the identity matrix and o$ is the z-component 
of Pauli matrix on C 2 . To find a proper value to we 
finally measure the value of 03 in the state p m such that 

M m ee trp m a 3 . 

We obtain 
Theorem 9 

Pm = {I + fm 2 iq2) ° 3 \ an d M m = .r(q 2 ). 



Thus the question is whether we can find such a to in 
polynomial steps of n satisfying the inequality M m > i 
for very small but non-zero q 2 . Here we have to remark 
that if one has q — then p — P and we obtain M m = 
for all m. If q ^ 0, the stochastic dynamics leads to the 
amplification of the small magnitude q in such a way that 
it can be detected as is explained below. The transition 
from po to p m is nonlinear and can be considered as a 
classical evolution because our algebra generated by Po 
and Pi is abelian. The amplification can be done within 
at most 2n steps due to the following propositions. Since 
f m (q 2 ) is x m of the logistic map x m+ i = f(x m ) with 
x Q = q 2 , we use the notation x m in the logistic map for 
simplicity. 

Theorem 10 For the logistic map x n+ \ = ax n (1 — x n ) 
with a € [0, 4] and x$ € [0, 1] , let xo be and a set J be 
{0, 1, 2, • • • , n, • • • 2n} . If a is 3.71, then there exists an 
integer m in J satisfying x m > i . 

Theorem 11 Let a and n be the same in the above 
proposition. If there exists too in J such that x mo > \ , 
then to > j^jj- 

According to these theorems, it is enough to check 
the value x m (M m ) around the above too when q is 
for a large n. More generally, when q—-^: with some 
integer k, it is similarly checked that the value x m (M m ) 
becomes over ^ within at most 2n steps. 

The complexity of the quantum algorithm for the SAT 
problem was discussed in Section 3 to be in polynomial 
time. We have only to consider the number of steps in 
the classical algorithm for the logistic map performed on 
quantum computer. It is the probabilistic part of the 
construction and one has to repeat computations several 
times to be able to distinguish the cases q = and q > 0. 
Thus it seems that the quantum chaos computer can solve 
the SAT problem in polynomial time. 

In conclusion of [10], the quantum chaos computer 
combines the ordinary quantum computer with quantum 
chaotic dynamics amplifier. It may go beyond the usual 
quantum Turing algorithm, but such a device can be pow- 
erful enough to solve the NP-complete problems in the 
polynomial time. The detail estimation of the complexity 
of the SAT algorithm is discussed in [21]. 

In the next two sections we will discuss the SAT prob- 
lem in a different view, that is, we will show that the same 
amplification is possible by unitary dynamics defined in 
the stochastic limit. 

IV. QUANTUM ADAPTIVE SYSTEMS 

In this section we begin to develop our programmeof 
constructing a physically realizable quantum amplifier for 
the OM quantum-SAT algorithm which is entirely within 
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the frame of standard quantum computation, namely: 
unitary evolutions implemented by the usual physical in- 
teractions (in fact we will consider the simplest class of 
these interactions: the dipole type ones). 

A classical amplifier is a sensor which reacts differently 
to the state of the input system, in other words it is an 
adaptive system. In our case the input state is the output 
of the OM algorithm, which is a quantum superposition. 
Thus we need a quantum amplifier and this naturally 
leads to the problem of developing a physically ealizable 
notion of quantum adaptive system. 

The idea to develop a mathematical approach to adap- 
tive systems, i.e. those systems whose properties are in 
part determined as responses to an environment [6,24], 
were born in connection with some problems of quantum 
measurement theory and chaos dynamics. 

The mathematical definition of adaptive system is in 
terms of observables, namely: an adaptive system is a 
composite system whose interaction depends on a fixed 
observable (typically in a measurement process, this ob- 
servable is the observable one wants to measure). Such 
systems may be called observable-adaptive. 

In the present paper we want to extend this point of 
view by introducing another natural class of adaptive sys- 
tems which, in a certain sense, is the dual to the above 
defined one, namely the class of state-adaptive systems. 
These are defined as follows: a state-adaptive system is a 
composite system whose interaction depends on the state 
of at least one of the sub-systems at the instant in which 
the interaction is switched on. 

Notice that both definitions make sense both for clas- 
sical and for quantum systems. Since in this paper we 
will be interested to an application of adaptive systems 
to quantum computation, we will discuss only quantum 
adaptive systems, but one should keep in mind that all 
the considerations below apply to classical systems as 
well. 

The difference between state-adaptive systems and 
nonlinear dynamical systems should be emphasized: 

(i) in nonlinear dynamical systems (such as those 
whose evolution is described by the Boltzmann equa- 
tion, or nonlinear Schrodinger equation, . . ., ) the in- 
teraction Hamiltonian depends on the state at each time 
t: H T = H^pt) ■ Vt . 

(ii) in state-adaptive dynamical systems (such as those 
considered in the present paper) the interaction Hamil- 
tonian depends on the state only at each time t = 0: 
H I = H I ( Po ). 

The latter class of systems describes the following 
physical situation: at time t = — T (T > 0) a system S is 
prepared in a state -0-t and in the time interval [— T, 0] 
it evolves according to a fixed (free) dynamics £^[-t,o] so 
that its state at time is U[_Tfl]ip-T ='■ ipo At time t = 
an interaction with another system R is switched on and 
this interaction depends on the state ?Ao : H j = Hi(tp ). 

If we interpret the system R as environment, we can 



say that the above interaction describes the response of 
the environment to the state of the system S. 

Now from the general theory of stochastic limit [1] one 
knows that, under general ergodicity conditions, an in- 
teraction with an environment drives the system to a 
dynamical (but not necessarily thermodynamical) equi- 
librium state which depends on the initial state of the 
environment and on the interaction Hamiltonian. 

Therefore, if one is able to realize experimentally these 
state dependent Hamiltonians, one would be able to drive 
the system S to a pre-assigned dynamical equilibrium 
state depending on the input state tpo- 

In the following section we will substantiate the general 
scheme described above with an application to the quan- 
tum computer approach to the SAT problem described 
in the previous sections. 

V. STOCHASTIC LIMIT AND SAT PROBLEM 

We illustrate the general scheme described in the pre- 
vious section in the simplest case when the state space 
of the system is Tts = C 2 . We fix an orthonormal basis 
of Hs as {e , ei}. 

The unknown state (vector) of the system at time t = 

ip ■= ^2 a e e £ = a oe + «iei ; ||-0|| = 1. 

££{0,1} 

In the case of Sec. 3, aicorresponds to q and ej does to 
Ii) (j = 0, 1) . This vector is taken as input and defines 
the interaction Hamiltonian in an external field 

H I = \\1>)y,\®(A+ + A g ) 

= Xa £ a £ \e £ }(e £/ \ ® (A+ + A g ) 

where A is a small coupling constant. Here and in the 
following summation over repeated indices is understood. 

The free system Hamiltonian is taken to be diagonal 
in the e e -basis 

H s ~ E s \e £ )(e s \=E \e )(e \+E 1 \e 1 )(e 1 \ 

££{0,1} 

and the energy levels are ordered so that E < E\. Thus 
there is a single Bohr frequency uj ■— E\ — E > 0. The 
1-particle field Hamiltonian is 

S t g(k) = e^g(k) 

where ui{k) is a function satisfying the basic analytical as- 
sumption of the stochastic limit. Its second quantization 
is the free field evolution 

e itH 0A ± ge -itH = A ±^ 

We can distinguish two cases as below, whose cases cor- 
respond to two cases of Sec. 3, i.e., q > and q = 0. 



G 



Case (1). If ao,ai ^ , then, according to the gen- 
eral theory of stochastic limit (i.e., t — > t/A 2 ) [1], the 
interaction Hamiltonian Hi is in the same universality 
class as 

Hi = D® A+ + D + ® A g 

where D := |e )(ei| (this means that the two interac- 
tions have the same stochastic limit). The interaction 
Hamiltonian at time t is then 



Hi(t) = e- itu>a D®A + s + h.c. = D®A+(e 



it(u(p)—u!o) 



g)+ h.c. 



and the white noise ({6*}) Hamiltonian equation associ- 
ated, via the stochastic golden rule, to this interaction 
Hamiltonian is 

d t U t = i(Db+ + D+b t )U t 

Its causally normal ordered form is equivalent to the 
stochastic differential equation 

dU t = {iDdB+ + iD + dB t - 1 _D + Ddt)U t , 

where dB t := btdt. 

The causally ordered inner Langevin equation is 

dj t (x) = dU* t xU t + U^xdUt + dU* t xdU t 

= Ut(-iD+xdB t - iDxdB+ - ~i_D+Dxdt + ixDdB^ 

+ixD+dB t - j-xD+Ddt + j-D+xDdt)U t 
= i Jt ([x, D+])dB t + ij t ([x, D])dB+ 

-(Re 7 _)j t ({P+P, x})dt + i{Imj-)j t ([D+D, x])dt 

+j t {D+xD){Re 1 _)dt, 

where j t (x) := UfxU t - Therefore the master equation is 

j t P\x) = {Im^)i[D+D,P\x)] - (Re-f-){D+D,P t (x)} 
+ (Re 7 _)P>+P*(x)P> 

where D + D = |ei)(ei| and D + xD — (eo, xe )|ei)(ei|. 

The dual Markovian evolution P' acts on density ma- 
trices and its generator is 

L,p = (Im^)i[p,D+D}-(Rc 7 _){p, D+D}+(Rc j_)DpD+ 
Thus, if po = |eo)(eo| one has 

L*po = 

so po is an invariant measure. From the Fagnola- 
Rebolledo criteria (cf. [25]), it is the unique invariant 
measure and the semigroup exp(tL^) converges exponen- 
tially to it. 

Case (2). If a\ = 0, then the interaction Hamiltonian 
Hi is 

Hi = A|e )(e | <g> (A+ + A g ) 



and, according to the general theory of stochastic limit, 
the reduced evolution has no damping and corresponds 
to the pure Hamiltonian 

H s + |e )<e | = (P + l)|e ><e | + Pi|ei)(ei| 

therefore, if we choose the eigenvalues E\,Eq to be in- 
tegers (in appropriate units), then the evolution will be 
periodic. 

Since the eigenvalues Pi, Po can be chosen a priori, by 
fixing the system Hamiltonian Hs, it follows that the pe- 
riod of the evolution can be known a priori. This gives 
a simple criterium for the solvability of the SAT prob- 
lem because, by waiting a sufficiently long time one can 
experimentally detect the difference between a damping 
and an oscillating behavior. 

A precise estimate of this time can be achieved either 
by theoretical methods or by computer simulation. Both 
methods will be analyzed in the full paper [3] . 



VI. CONCLUSION 

We showed in [7,9,10] that we could find an algorithm 
solving the SAT problems in polynomial steps by com- 
bining a quantum algorithm with a chaos dynamics. We 
used the logistic map there, however it is possible to use 
other chaotic maps if they can amplify one of two coffi- 
cients. In this short paper we pointed out that it is possi- 
ble to distinguish two different states, yl — q 2 |0) +q |1) 
(q 7^ 0) and |0) by means of the adaptive dynamics and 
the stochastic limit. Finally we remark that our algo- 
rithm can be described by deterministic general quan- 
tum Turing machine [22,23], whose result is based on the 
general quantum algorithm mentioned in Sec. 2. 
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